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An electronic phase separation model provides a natu- 
ral explanation for a large variety of experimental results in 
the cuprates, including evidence for both stripes and larger 
domains, and a termination of the phase separation in the 
slightly overdoped regime, when the average hole density 
equals that on the charged stripes. Several models are pre- 
sented for charged stripes, showing how density waves, super- 
conductivity, and strong correlations compete with quantum 
size effects (QSEs) in narrow stripes. The energy bands asso- 
ciated with the charged stripes develop in the middle of the 
Mott gap, and the splitting of these bands can be understood 
by considering the QSE on a single ladder. 



I. INTRODUCTION 

Possible electronic phase separation (EPS) in the 
cuprates has usually been found in terms of stripe phases. 
Thus, neutron diffraction measurements find evidence for 
fluctuating stripe order in La2-xSra;Cu04 (LSCO) associ- 
ated with incommensurate inelastic neutron scatteringEI, 
whickjCan be transformed into long range charged stripe 
ordero by co-doping with Nd or Eu. Similar incommen- 
surate peaks are found in other cuprates, and long range 
charge ordered states are also foundd in strongly under- 
doped samples of YBa2Cu307_5 (YBCO), with stripes 
parallel to the chains, while at higher doping short range 
stripe order is found at virtually all temperatures up to 
the pseudogap T*. However, EPS can also manifest itself 
in the form of larger domains, particularly if the dopant 
ions are somewhat mobile and can follow the hole mo- 
tion. These domains, long known in La2Cu04+5 may re- 
cently have been found in. scanning tunneling microscope 
(STM)cl and microwaveS studies of Bi2Sr2CaCu208+5 
(Bi2212). 

Stripes are typically interpreted either in terms of 
charged domain walls in antiferromagnets or of restricted 
EPS. While the models are virtually indistinguishible 
in the underdoped regime, they diverge as doping in- 
creases. A phase separation model implies the existence 
of a uniform end phase: at some doping antiferromag- 
netic (AFM) stripes must disappear, leaving a uniform 
state similar to that found on the charged stripes. For 
slightly lower doping the AFM stripes form a minority 
- a situation not possible in a domain wall picture. (A 
note of caution: more recent domain wall models find 
a lower doping on the charged stripes, and as the dop- 
ing decreases, the distinction between the two models 
gradually blurs.) Evidence fqs the termination of the 
AFM stripes at a fixed dopujigEl and recent observations 
of larger domain structuresQ taken together provide ex- 



tremely strong support for a phase separation scenario. 

At this point it is essential to better characterize the 
nature of the stripes, in particular the charged stripes, 
and to understand how their properties affect pi. 
cal phenomena, injjjarticular photoemission spectra.! 
While fluctuationsElEl play an important role in the real 
cuprates, we have constructed ordered stripe arrays, for 
which detailed tight-binding calculations are possibleO. 
Such calculations can aid in elucidating the structure of 
the charged stripes, both for wide stripes ~ what stabi- 
lizes the preferred hole density - and for narrow stripes 
- how do quantum size effects (QSE) modify properties 
of the stripes. The present paper provides an extended 
analysis of thesejssues; some of these results were sum- 
marized recentlyli3. 

The paper is organized as follows. Section II enumer- 
ates key issues which must be addressed in any EPS 
model of stripes, including determination of the hole 
density on charged stripes. It is found that the model 
can simultaneously account for a large number of exper- 
imental results. The charged stripes are probably stabi- 
lized by some competing order, either magnetic or para- 
magnetic (including charge density waves). Section III 
shows that magnetic charged stripes can arise in a mean 
field Hubbard model, and can be either ferromagnetic 
or a linear antiferromagpstic (LAF) phase similar to the 
White-Scalapino stripescj. Calculations on stripe arrays 
find that the charged stripes lead to midgap states near 
the Fermi level; Section IV shows how these data can 
be interpreted in terms of quantum size effect (QSE) 
on single stripes. Section V presents the results of sin- 
gle stripe calculations: competing charge density wave 
(CDW) and superconducting orders can exist on para- 
magnetic stripes, but they are strongly modified by the 
QSE. On LAF stripes, d-wave superconductivity and an 
unusual form of CDW are both found to persist down to 
the narrowest (2 cells wide) stripes. Section VI includes 
results on stripe arrays: the LAF stripes produce pho- 
toemission constant-energy maps in significantly better 
agreement with experiment. In the Discussion, Section 
VII, we summarize additional recent evidence which fa- 
vors a phase separation model - in particular, evidence 
that AFM stripes terminate slightly beyond optimal dop- 
ing, where superconducting properties remain strong, 
and that a second regime of phase separation exists above 
optimal doping - and we show that experiments are con- 
sistent with our prediction that the superconducting gap 
grows in the underdoped regime. Conclusions are given 
in Section VIII. 
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II. KEY ISSUES FOR AN EPS MODEL OF 
STRIPES 

A phase separation model of stripes is characterized 
by the two well-defined stable end phases between which 
phase separation takes place. The insulating stripes are 
generally understood to be antiferromagnetic (AFM) - 
essentially the same as the Mott insulator found in un- 
doped cuprates. The hole-doped stripes are assumed to 
have a finite doping, xq ■ This simple idea has three exper- 
imentally verifiable characteristic features: (i) the stripe 
phase must terminate when x — xq] (ii) there will be a 
crossover at a lower doping, Xcr ~ xq/2 from magnetic- 
dominated {x < Xcr) to charge-dominated (x > Xcr) 
stripe arrays; (iii) some interaction on the charged stripes 
stabilizes the end phase at xq. Here we discuss current 
evidence for (a) the doping on the charged stripes, (b) ev- 
idence for a crossover, and (c) the nature of the dominant 
interaction stabilizing the charged stripes. In addition, 
we ask what constraints the Yamada plot puts on the 
model, and how a model based on EPS compares to a 
domain wall model. 



A. Hole Doping on Charged Stripes 

Recent evidence suggests that the stripes and pseudo- 
gap terminate at the same dopiiie; xo while superconduc- 
tivity persists to higher dopinga. However, the proper 
choice of Xq requires some discussionj-. The neutrop. 
diffraction measurements of Tranquadaa and YamadaEJ 
have established that charged stripes in La2-2:Sr^Cu04 
(LSCO) have an invariant topology over the doping range 
0.06 < X < 0.125, acting as antiphase boundaries (APBs) 
for the AFM stripes and having a net doping of 0.5 holes 
per stripe. However, there are two models for how this 
charge is distributed: either in one row with average hole 
density 0.5 hole per copper site or in two rows with 0.25 
hole per copper. These two alternatives are often some- 
what simplistically referred to as site order vs bond or- 
der. The strongest evidence distinguishing between tha 
alternatives comes from x-ray data on the charge ordertil 
at x=l/8: non-observation of diffraction harmonics sug- 
gests a sinusoidal distribution of charge. For a 4-Cu 
repeat distance, two insulating and two charged rows 
would be exactly sinusoidal, whereas one charged and 
three insulating rows should have signifigant harmonic 
content. A similar conclusion was reached by /kSR line- 
shape analysisEj. However, the charge ordering peaks are 
weak, and it remains possible that fluctuations or disor- 
der could wash out the harmonics. 

Direct evidence for the density on the charged stripes 
is found from low temperature NQR measurementsli^, 
which find values xq ^ 0.18 — 0.19. While this is close to 
the lower value, the small difference can also be under- 
stood: this is a ZocaZ measurement, and it is expected that 
some holes will be pushed off onto the magnetic stripes. 



Indirect evidence favoring the lower hole density includes 
the fact that it is easier to understand the properties of 
AFM stripes in terms of even-leg ladders (e.g-i the AFM 
stripes would be two coppers wide at x=l/8)EZl, and that 
the stripe phase appears to terminate when the average 
doping approaches a;=0.25 (Section VII.A.l). The lower 
doping is also more consistenj-jvith the tJ model simula- 
tions of White and ScalapincO. 

Tallonlla finds an optimal doping at Xopt = 0.16 for 
all cuprates, with respect to which the stripes terminate 
at a doping x = 0.19. However, it is hard to recon- 
cile a^common optimal doping with muon spin resonance 
datallj, which find Tc is optimized at very different val- 
ues of Us/m iris is the paiii-density, which seems to scale 
with carrier concentrationcj, and m effective mass) for 
LSCO and YBCO. We assume instead that Xopt scales 
with n/m, so if Xopt — 0.16 for LSCO, it is |£L21 for 
YBCO, in good agreement with several cstimatcfH. This 
also resolves a problem with the thermopower. While the 
thermopower appears to be universal for most cuprates, 
and the best means of estimating the doping is from the 
room temperature thermopower, LSCO is anomalouSrin 
that 'overdoped' samples still have high thermopowerc3. 
If the doping for YBCO is rescaled as above, however, the 
thermopower data of LSCO fall on the universal curve. 
Hence, the anomaly for LSCO is not in the thermopower, 
but in a too low value of Tc, which is accompanied by a 
too low value of Xopt ■ It is likely that these features are 
associated with a competing LTT phase, which is most 
prominent in LSCO, and which also leads to the most 
nearly static stripe correlations. 

Taking a;opt=0.21 for YBCO gives xq = 19/16 x 0.21 = 
0.25, which we believe holds for all cuprates, includ- 
ing LSCO, Section VII.A.l. This would lead to very 
wide charged stripes near optimal doping: the width of 
the charged stripes N satisfies N/{N + 2) = 16/19, or 

= 32/3 ~ 10 Cu wide. Hence, models of isolated quasi- 
one-diniensional charged stripes are likely to be valid only 
in the far underdoped regime, while for the good super- 
conductors^ better model would be a metal with intrinsic 
weak linkscj. 



B. Crossover at 1/8 Doping 

For Xq ~ 1/4, the crossover Xcr = 2;o/2 can be identi- 
fied with the 1 /8 anomaly, where both charged and AFM 
stripes have their minimal width (2 Cu atoms). There is 
considerable evidence that the doping 1/8, in addition to 
its special stability, acts as a crosscper in the properties 
of the stripes. Thus, Uchida, et al.Eil, studying the Hall 
coefficient Rh^ find a crossover from one-dimensional be- 
havior {Rh — > as T 0) for x < Xcr = 1/8 to two di- 
mensional behavior (coupled charged stripes) for x > Xcr- 
In YBCO, the spin gap grows slowly with doping for 
X < Xcn then more rapidly for x > Xcr] this behavior 
can be understood in terms of coupled spin ladders, as the 
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coupling changes with the, width of the charged stripesB. 
In Eu substituted LSC0l3, the Meissner fraction is neg- 
ligibly small for x < Xcr, then grows roughly linearly 
with doping until x ~ 0.18, staying large up to at least 
X = 0.22. Finally, the two-magnon Raman peak in LSCO 
has a splitting at |lo|W temperatures which has been associ- 
ated with stripesc3, on analogy with similar observations 
in La2-a:Sra:Ni04 (LSNO)EI. For x < Xcr the ratio of 
the two peak frequencies is constant and consistent with 
a simple stripe model; for x > Xcr the lower frequency 
starts decreasing with doping. Moreover, the higher fre- 
quency loses intensity with doping; near x — 0.26, the 
intensity of one mode approximately disappears, while 
the frequency of the other mode extrapolates to zero. 



C. Constraints from the Yamada Plot 
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FIG. 1. (a) Yamada plot of incommensurability 5 vs doping 
X for LSCO. Open ssuares = elastic neutron scattering in Nd 
substituted samplesH; others =ji inelastic neujtspn scattering 
for vertical stripes (opep-icirclesa or diamondfj) or diagonal 
stripes (filled diamondstj). Dashed (solid) line = prediction 
of EPS model without (with) commensurability effect at 1 /8 
doping. (b,c) = Stripe phase model without commensurabil- 
ity effect, at a; = 1/12 (b) and 1/6 (c). For this figure, the 
charged stripes are assumed to have linear antiferromagnetic 
order (Section III). (d,e) = domains associated with commen- 
surability pinning of 1/8 doped phase; (e) = blowup of (d). 

The YamadaElSS plot. Fig. |l|a, provides a severe con- 
straint on any model of stripes: in LSCO the incom- 
mensurability b is found to grow linearly with doping for 
X < Xcr but to saturate for x > x^r- Furthermore, the 
saturation value is just the incommensurability expected 
for X = I/8-doped stripes, 5sat = Xcr- A similar satu- 
ration has been reported in YBCO, but different groups 
find different values for 5sat- ^ l/6[ |9| or - I/IO[ |o). 



If one of these values proves correct, it would suggest 
some nonuniversality in xq, perhaps associated with bi- 
layer splitting. In the domain wall model, the saturation 
is interpreted as evidence that the stripes stop changing 
in width, and additional holes leak into the antiferromag- 
netic background, washing out the stripes. 

For X < 1/8, the EPS model agrees with the Yamada 
plot. Fig. I^b: increase in doping causes the AFM stripes 
to narrow, with no changes in the charged stripes. Note 
that, for concreteness we have assumed that the charged 
stripes have LAF order (Section III); these stripes natu- 
rally act as antiphase boundaries (APBs), consistent with 
the neutron evidence. However, the naive prediction of 
the EPS model is in disagreement with the Yamada plot 
for X > 1/8, Fig. 0c: as the charge stripes get wider, the 
incommensurability should decrease, while the neutron 
peaks may broaden if the wider charge stripes do not act 
as APBs. This behavior is not observed. However, the 
model can be simply modified to explain the observed 
saturation. Fig. ^,e. This would be a commensurabil- 
ity effect, with part of the sample pinned at 1/8 doping 
while the rest forms a different phase (e.g., at 1/4 dop- 
ing) where no stripes are present. Such behavior is well 
known in nickelates, where coexistence of 1/3 and 1/2 
stripes is common. 

Evidence for such pCpmmensurability effects can be 
found _by comparipgta the chemical potential jj, in 
LSCOEl and LSNOB. In both materials, /i IS approxi- 
mately independent of doping between half filling (x — 0) 
and Xcr, with Xcr = 0.125 in LSCO, 1/3 in LSNO, which 
has long-range charge order. Remarkably, the simple pre- 
diction of a macroscopic phase separation model (/z con- 
stant throughout the phase separation regime) is violated 
in LSNO: fj, is not constant for 1/3 < a: < 1/2, although 
the stripe phase persists over the full doping range. This 
is presumably a commensurability effect, reflecting the 
coexistence of 1/3 and 1/2 stripe phases for intermedi- 
ate dopings. The same doping dependence of fi is found 
in LSCO: constant for a; < 1/8, variable for x > 1/8. 
Indeed, the LSNO data can be scaled to that of the 
cuprates, giving a cammon doping dependence fj,{x/xc) 
in both compoundgll3. 



D. Domain Phase 

Recently local charged domains have been directly vi- 
sualized by STM studies in Bi2212a with gaps ranging 
from slightly overdoped to strongly underdoped. Such 
behavior is extremely difficult to explain in the domain 
wall picture, but is easy to understand in terms of EPS, 
and indeed is consistent with the commensurability ef- 
fects discussed above. The gap distribution is found to 
be broad rather than bimodal, but that is expected since 
EPS is very sensitive to charged impurities, and the local 
gap will correlate with the local density of dopant, pre- 
sumably interstitial oxygen in Bi2212. This sensitivity 
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to impurities leads to the question of which came first: 
is a preexisting EES sensitive to local impurities, or does 
oxygen clusteringE3 provide the driving mechanism for 
domain formation? Since the electronic inhomogeneity 
seems characteristic of most cuprates while there is con- 
siderable variety in the doping counterions, the simpler 
interpretation would appear to be that the EPS is pri- 
mary. Thus, in La2Cu04+5 (LCO), the interstitial oxy- 
gens are highly mobile, allowing the domains to grow to 
macroscopic size. Similar clusters form in YBCO (here 
associated with chain oxygens), but carv-he suppressed in 
fully oxygenated samples by quenchingS^. On the other 
hand, well formed stripes appear when the doping coun- 
terions are least mobile, in LSCO. 

The domains in Bi2212 would seem to be consistent, 
but a number of questions remain about the role of an- 
nealing. While these domains are regularly found in STM 
studies when the samples are cleaved at low-temnerature, 
they can be annealed out at high temperaturesEEl. This 
could be caused by a 'melting' qLihe EPS, as is found 
in LCO near room temperature^j, and can be tested 
by careful annealing studies. Alternatively, it may be 
a question of pinning the EPS. Certainly, it is known 
that a weak domain disorder is necessar#-|to pin vortices, 
to observe the vortex lattice via STM.a Also, the role 
of the interstitial oxygen in the superlattice modulations 
in Bi2212 needs to be better understood. The strength 
of this modulation suggests that it is associated with 
an ordering of the interstitial oxygen, but this should 
lead to a strong correlation between the domains and 
the superlattice, which seems not to be the case. Fur- 
ther evidence for electronic inhomogeneity comes from 
microwave measurements: anomalies in Bi2212 ha^ie been 
interpreted in terms of similar electronic domainsQ, sug- 
gesting that they are representatiye of the bulk, while 
measurements on other cupratestD find similar anoma- 
lous behavior which was interpreted in terms of pinned 
CDW's, possibly stripe related. A domain picture would 
also explain the persistence of nodal quasiparticles,in the 
underdoped regime, at least down to 1/8 doping53. 



E. Comparison of Domain Wall Stripes and EPS 
Stripes 



Table I summarizes the discussion of the previous sub- 
sections, and compares the predictions of EPS models 

Table I: Comparison of Stripe Models 



with domain wall models of stripes. Domain wall 
models arose-, iftncarly unrestricted Hartree-Fock (UHF) 
calculationg53~Eil for the doped tJ model, as heteroge- 
neous ground states in which the holes are confined to 
domain walls between AFM domains which act as APBa. 
These domain walls are not driven by phase separationcj; 
the phase is realized as a long-period modulated AFM, 
with holes doping the rows of spins where the moment 
changes signE3. When t' = 0, there is one hole per row, a 
large hole doping which is inconsistent both with experi- 
ments on stripes in the cuprates, and with depsity matrix 
rcnormalization group (DMRG) calculationst2l. More re- 
cent calculations reduce the doping, in ways that bring 
the model close to the phase separation picture: (1) as- 
suming CDW order along the charged stripea^J lowers xq 
to 0.5, while (2) letting t' < results in xq ~ 0.2.N Since 
the CDW order (case 1) lowers the free energy of the 
charged stripes, it should be possible to dope the system 
all the way to a pure CDW at x = 0.5. In case 2, it is not 
explained why a low density appears; a good possibility 
is that this doping corresponds to the Van Hove singu- 
larity (VHS), as in the phase separated stripes. (This 
could be checked by varying t' and calculating X[){t').) 
If such low-density domain walls do exist, the striped 
phase must somehow terminate when the average doping 
approached 0.2, but howrtiiis happens is not explained. 
However, this calculationcJ is not fully UHF, bemg re- 
stricted to periodic arrays, and UHF calculationsEl find 
instead ferromagnetic charged stripes which can be un- 
derstood within a phase separation model (see below). 

Since the stability of these t' 7^ domain wall stripes 
is unclear, in Table I a comparison is made with the ear- 
lier, higher-hole-density domain wall models. For such 
high doping, stripe termination at x = 1/4, or a major 
crossover near 1/8, are both difficult to interpret. Stripes 
as APBs are more natural in domain wall models, but we 
shall show that they can also arise in EPS models. Note 
that the distinction between the models may be becom- 
ing moot, since in the domain wall rnodels there is now a 
search for a 'hidden' order parameterca, and it is common 
for secondary order parameters to be expressed predom- 
inantly on domain walls of the primary orderEl. 

In conclusion, assumption of a charged stripe doping 
xo — 0.25 reconciles the neutron diffraction data, evi- 
dence for a termination of the stripe phase near xq, and 
the 1/8 anomaly as a crossover effect near while 
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commensurability effects can explain tlie saturation in 
the Yamada plot and the STM observation of charge do- 
mains. Further supportive evidence will be presented in 
Section VILA. In the remainder of this paper, we will 
assume that xq = 0.25. 



F. What Stabilizes Charged Stripes? 

In any phase separation model, a key issue is under- 
standing the nature of the charged stripes. Indeed, since 
superconductivity seems to arise predominantly on these 
stripes, such understanding is likely to play an important 
role in elucidating the origin of the high superconducting 
transition temperatures. For the stripe phase to exist, 
the doping a;o must be particularly stable. This can arise 
via an electronic instability, which opens up a gap over 
much of the Fermi surface, making the electronic phase 
nearly incompressible. This 'Stability from Instability' is 
a fair bi_ general feature, underlying, e.g., Hume-Rothery 
alloysB JThis is a modification of an argument due to 
Andersoric3.] Here, we explore a number of candidates 
for the predominant jsLectronic instability. 

In a related paperE^, we will provide strong evidence 
that this 'hidden order' is a form of CDW, which could 
include the flux phase. However, here we will explore 
a wider variety of possibilities. One issue is that in 
the low doping limit the charge stripes act as APBs for 
the AFM stripes. Such an effect arises naturally if the 
charged stripes have some residual magnetic interaction, 
and we will explore this possibility. However, in nick- 
elates charged stripes coupled to a CDW are found to 
act as APBalj. The large Hubbard on-site repulsion U 
plays an important role. Strong correlation effects lead 
to two classes of charged stripe phases: either the con- 
straint against double occupancy leads to magnetic or- 
der (magnetic charged stripes) or kinetic energy domi- 
nates, leading to a magnetically disordered phase (para- 
magnetic charged stripes) j— with correlations Jeading to 
reduced hopping, as in tJcJ and slave bosorca models. 
Section HI will provide examples of both classes, denoted 
as Class B and Class A stripes, respectively. Class A 
stripes could be simply correlated metals (as in tJ or 
slave boson calculations) or could have additional, e.g., 
CDW, order. A crossover from magnetic to correlated 
paramagnetic groundstate arises aSpS, function of doping 
in models of itinerant ferromagnetsca. 

For completeness. Class C stripes are defined as those 
arising not from phase separation but from long-range 
modulated AFM order (domain wall stripes). Such 
stripes have been described in detailE3, and will not be 
further considered here. 

In the next Section, we show that Class B (magnetiii) 
charge stripes can arise in a mean-field Hubbard modeO, 
with the charged stripes displaying either ferromagnetic 
(FM) or linear antiferromagnetic (LAF) (ordering vector 
(tt, 0)) order. The LAF stripes are very similar to White- 



Scalapino stripeslS. The FM phase is stabilized by VHS 
nesting. This FM phase may be present in ruthenatesEil, 
but is unlikely to be relevant for the cuprates (for one 
thing, the FM^stripes are likely to be diagonal, and do 
not form APBscil, contrary to experiment). We have sufe 
gested that other VHS-stabilized phases are more likelyEj 
(see also ReL 53), and here we explore the properties of 
a Class A CDW phas(£3. At a doping x ~ 0.25, the 
effects of correlations are relatively weak, renormalizing 
the bandwidth by a factor of ^2. Thus for the present 
calculations on paramagnetic stripes, it will be assumed 
that renormalized parameters are used, and other effects 
of strong correlations will be neglected. 

While LAF stripes are most stable when t' = 0, we 
explore the possibility that they can be stabilized even 
when i' 7^ by on-stripe CDW or superconducting order. 
Ordinarily, the CDW phase is believed to compete with 
strong correlation effects, but we find (Section V.B) that 
an unusual form of CDW phase can coexist with LAF 
order: the charge density varies between zero and one 
(not two) holes per atom. 



G. Notation on Stripes 

In a stripe array, the alternating stripes are associ- 
ated with the two stable thermodynamic phases. Here, 
we summarize the different ways these stripes are de- 
noted in this paper. The stripe with lower hole doping 
is variously referred to as 'insulating' or 'antiferromag- 
netic' (AFM). (These are also known as magnetic stripes, 
since the charged stripes have considerably weaker mag- 
netic order, but we will avoid that notation here.) The 
stripe with higher hole doping is generically referred to 
as 'charged' or 'hole-doped'. At low temperatures, these 
stripes are also 'superconducting', but at high temper- 
atures, they are stabilized by some 'hidden order', and 
one purpose of this paper is to explore a number of pos- 
sible orders. The orders fall into two groups: 'magnetic 
charged' (Class B) stripes could have FM or LAF or- 
der (the latter are White-Scalapino-like stripes), while 
'paramagnetic' (Class A) stripes could have CDW or flux- 
phase order. 



III. PHASE SEPARATION IN A MEAN-FIELD 
HUBBARD MODEL 

Strong coupling models would seem to be natural for 
producing phase separated or striped ground states. Any 
magnetic ordering avoids double occupancy, while chang- 
ing from one form of magnetic order to another, via, e.g., 
doping, requires highly collective spin rotations, as com- 
peting orders are orthogonal. While superexchange leads 
to antiferromagnetic (AFM) insulators at half filling, 
doping tends to favor textures with parallel spins (e.g., 
ferromagnetic (FM)) to maximize hole hopping. Such 
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ferron phases were introduced long before high TcEJ, but 
it remains controversial whetiier such states are ground 
states of the Hubbard modelEa. While the t J model does 
have phase separation for large J/t, it is unclear whether 
such phases |extciyLto the values J jt ^ 0.3 expected for 
the cupratesBETEl. 

We have found phase separated solutions of the Hub- 
bard model at mean-field leveEJ. While these solutions 
are metastable in UHF calculations, they closely resem- 
ble the WS stripes, and provide an interesting example of 
phase-separation mediated stripe phases. We find a well 
defined surface tension for wide, isolated stripes, which 
decreases and changes sign as the stripes become nar- 
rower. When the surface tension becomes negative, the 
stripes no longer remain straight, and spontaneously me- 
andering solutions are found. 

These solutions are found by considering only low- 
order commensurate phases, with wave vector qx,qy ^ 
01 Qi — n /a only. The bare dispersion is 



Efc = -2t{Cx + Cy) ~ At'CxC 



1" 



(1) 



with Ci = cos/cifl. The Hubbard interaction C/^j(rti| — 
1/2) (nil — 1/2) leads to magnetic order with a mean- field 
magnetization vriq at wave vector and the quasiparticle 
dispersion becomes 



±Eo, 



where 



and 



(2) 



(3) 




FIG. 2. Free energy vs. doping for several magnetic phases 
of the Hubbard model assuming U — 6.03t, and t' = — 0.276t. 
Diamonds — AFM, triangles — LAF, circles = FM, and 
squares — PM phase. Dashed lines = tangent construction; 
dot-dashed line — Eq. ^. Inset: Dispersion of FM phase 
at a; = 0.31; Brillouin zone points F = (0,0), X = (7r,0), 
S = (7r,7r). 

The competing phases include AFM ioi q — Q = 
(tt, tt), FM with q — (0,0), and a linear antiferromagnet 
(LAF) with q = (tt, 0). When the LAF stripes are two 
cells wide, this LAF phase closely resembles the White- 
Scalapino stripes. Fig. |l|e. The AFM state has lowest 
free energy at half filling, but (for t' = 0) is unstable for 
finite hole doping. For t' = 0, there is phase separation 
between the AFM and LAF phases, while for finite t' the 
phase separation is between AFM and FM phases. Fig. ||. 
(When electron-phonon coupling is included, it is found 
that tije FM phase is unstable with respect to a CDW 
phaseES.) 



(4) 



For the cuprates, we expecto t ~ S25meV, U/t ~ 6 
and t' /t ~ —0.276. The magnetization is found self- 
consistently from 



m,=Y.^f{E^)-f{E+)) 



Uniq 



(5) 



with Fermi function f{E) = 1/(1 + e{E-Ep)/kBTy j.^^ 
resulting free energy is 



with 



F = Eg~TS + Uiml + —), 



fc,i=± 



(6) 



(7) 
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FIG. 3. Phase diagrams, x{U) for the Hubbard model, with 
t' = (a) or t' = -0.276f (b). (Dashed line in (b) = doping 
of VHS.) 



S = kBj2 U{E{)lnU{E,)) + (1 - f{E,))ln{l - f{E,))). 

(8) 



k,i=± 



The LAF stripes for i' = are discussed in Ref . For 
t' ^ electron-hole symmetry is absent; for hole doping 
x > 0, there is phase separation to a FM phase, con- 
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sistent with recent simulations by Vozmediano, et al.0, 
who find a uniform FM phase at x = 0.15 for U = 8t, 
t' = 0.3t. However, on the electron- doped side a uniform 
AFM phase is stable over a large doping range, suggestive 
of the asymmetry found in the cuprates. The dot-dashed 
curve shows that in the phase separation regime, the low 
energy physics can be approximated by the form of free 
energy assumed in Ref. 0, 



fix) = /ox(l 



(9) 



(neglecting a term linear in x) with xq the hole doping of 
the uniform charged phase. The FM phase is stabilized 
by VHS nesting, inset to Fig. ^, as found previouslytil. 
The tangent construction tends to select the FM phase 
at dopings somewhat away from optimal nesting (Fermi 
energy above midgap). Both regimes of phase separation 
seem to be driven by hole delocalization: one-dimensional 
(along the LAF rows) when t' = 0, two-dimensional for 
finite t' . 

The resulting phase diagrams x ys U , Fig. ^, are strik- 
ingly different. For t' = 0, Fig. the phase separation 
is between the AFM and paramagnetic (LAF) phase for 
U < Uc = 5.3i {U > Uc), while for finite t' , Fig. |b, 
there is generally a VHS-stabilized FM phase. For small 
U and t' 7^ 0, there is a regime where simple spin-density 
wave theory works and a uniform AFM phase is stable, 
but when t' — phase separation persists for all finite 
U. Note that Uc marks a crossover between Class A 
(paramagnetic) and Class B (magnetic) charged stripes. 
The value Uc is close to the U = 6.03t expected in the 
cuprates, although for finite t' Uc decreases, Uc ~ 3t for 
t' = — 0.276f, and the range of U for which paramag- 
netic stripes are stable becomes very small. While these 
stripes are metastable in UHF calculations, we will show 
below that the stability of phase separating stripes can 
be enhanced by additional interactions beyond the pure 
Hubbard model. 



IV. ISOLATED STRIPES VS. ARRAYS 

In ordered stripe arrays it is foundi that the AFM 
stripes have a Mott-Hubbard gap, and the features near 
the Fermi level are associated with the charged stripes. 
This charge stripe dispersion shows a series of quasi-one- 
dimensional features which qualitatively resemble the 
bands of an isolated stripe produced by QSE. In this 
Section, we make a quantitative comparison with isolated 
stripes, and explore the mechanism of QSE-induced Van 
Hove splitting. 

For a single stripe N Cu atoms wide, the dispersion is 
still given by Eq. |l|, but the allowed values are quan- 
tized, with = km = mn / {N + 1), m = l,2,...,iV. 
These are in fact the usual quantized Bloch bands, but for 
large N the quantization is not noticeable. For small N, 



the dispersion appears as a series of overlapping one- 
dimensional (Id) dispersions. Equation |l| can be rewrit- 
ten as N Id dispersions 



2{t + 2t'c,n)Cy 



(10) 



These are the QSE, with corresponding density of states 
(dos) shown in Fig. 0. Notice that for N = 100, the 
VHS is readily apparent in the dispersion. Even down to 
N — 2, the VHS is clearly defined (albeit only within a 
finite interval) as the locus of energies where all subbands 
overlap. In fact, the QSE opens a gap at the VHS, ef- 
fectively lowering the kinetic energy of the electrons just 
like a conventional (e.g., CDW) gap. The VHS splitting 
can be found from Eq. 10 



VHS 



4t(l-ci). 



(11) 



This splitting has two consequences: first, the VHS split- 
ting enhances the stabilization of the striped phase; but 
second, the QSE gap competes with other gaps, such as 
CDW's and superconductivity. However, while the QSE 
splits the VHS peak, substantial dos remains ungapped, 
so additional instabilities remain possible. This compe- 
tition will be discussed further in the next Section. 
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FIG. 4. Density of states for a single stripe of width A'^ — 
100 (a), 10 (b), 4 (c), or 2 (d) atoms. Based on Eq. | with 
t7t=-0.276. 

It should be noted that this is the first direct demon- 
stration that the VHS can be defined on a stripe only 
two atoms wide. The definition is quite analogous to the 
standard definition in two dimensions: the point at which 
the bands cross over from electron-like to hole-like. 

Figure [s] compares the gaps of aijingle stripe with those 
found in the ordered stripe arrayO; the array is labelled 
(to, n) when the magnetic stripes are m coppers wide 
and the charged stripes are n coppers wide. In the array 
calculation, no competing order was introduced on the 
charged stripes, so the QSE provides the only gap. Fig- 
ure || shows that there is a very good match for both 2 Cu 
wide and 6 Cu wide stripes, although for the 2 Cu stripe. 



7 



the VHS gap is somewhat larger for the single stripe than 
in the array. From Eq. |l^, the Id bands have dos peaks 
at band bottom and top; the band bottom corresponds 
to ky = - i.e., the dispersion from F — > X is flat, at 
the energy corresponding to the lower dos peak. (The 
intensity along F — > X is given by a structure factor, 
which does not directly come into a single stripe calcu- 
lation.) Along X ^ S one should see the Id dispersion 
extrapolating to the band top at 5* = (tt, tt). Given the 
good agreement, it should be possible to analyze compet- 
ing orders on a single stripe, for which the calculations 
are simpler (no need to self-consistcntly adjust doping on 
each row to account for charging effects). 




FIG. 5. Dispersion of stripe arrays; (a) (6,2) array with 
charged stripes 2 Cu wide (AFM stripes 6 Cu wide); (b) (2,6) 
array with charged stripes 6 Cu wide (AFM stripes 2 Cu 
wide). Data from Figs. 7a and 7d respectively of Ref. ^ 
triangles (diamonds) = predominantly from charged (AFM) 
stripes, while circles — mixed origin; dashed line = Mott 
bands of AFM stripes; solid line — single (charged) stripe 
model, with approximated by nearest quantized value. 



V. ORDERING ON SINGLE CHARGED STRIPES 

A. CDW's and Superconductivity on Paramagnetic 
Stripes 

1. Electron- Phonon Coupling 

In this section, we will develop two closely related mod- 
els of the competition of CDW order and superconduc- 
tivity on a single charged stripe. The first is a Class 
A paramagnetic stripe, stabilized by electron-phonon 
couplingHo, while the second has dominant electron- 
electron coupling, with (Class B) magnetic charged stripe 
order. |— , 

The earlier calculation^ used a van Hove stabilized 
CDW model to describe the doping dependence of the 
pseudogap. Here we reapply the model for a single stripe, 
by introducing the following modifications. (1) The cal- 
culation is carried out on a single, finite width stripe. (2) 



Pinning to the VHS arises naturally, since all stripes are 
at the same doping. (3) For closer approximation to ex- 
periment, d-wave superconductivity is assumed. (4) Cor- 
relation effects due to Hubbard U are neglected: previous 
slave boson calculations suggest that the raain effect is a 
bandwidth renormalization of a factor ~2.cfl 

We briefly recall tha^ energy dispersion and the gap 
equations of the modelEj, generalized to d-wave. In terms 
of a function 



^ J 1, if - efI < TT-uJph] 
^ 1 0, otherwise, 



(12) 



the gap functions are A^: — /S.qQj^{cx — Cy)/2 for super- 
conductivity, and Gj: = Gq-\- Gi^k^^_^_Q for the CDW. 
The energy eigenvalues are E±^k and their negatives, with 



2Gl ± El), 



(13) 



El 



{El 



El+Q? 



^GM, E\ 



(Afc - Afe+g) , Cfei = £k± Sk+Q, and the nesting 



vector Q ~ (tTjTt). The gap equations are 



p2 ,2 
'^+,k-^k+Q 



A AaAEj 

e 



' El.k 
[Al + Gl] 



E^ 



E^ 



el - 

k+Q 



2E+,k 

'k+Q 



■ tanh 



X 

,k 

E+,k 



2E. 



tanh 



2kBT 

E-.k 



,k 



2kBT 



(14) 



E+M + '^k'^k+Q 



G^k^ 



e.Gj: 



E- 



2E- 

E-,k + ^k^k+Q 



— tanh 



-,k 

E+^k 



+,fc 



2E^ 



- Gl 
— ^tanh 



2kBT 

E-h ' 



,k 



2kBT 



(15) 



with interaction energies Aa and Ag, and Oq — ^k^k+Q^ 
Gl = 1. 

The CDW-superconducting competition was studied 
in bulk in Ref. |6^. The previous results are recovered for 
a sufficiently wide stripe (~100 Cu wide). For narrower 
stripes, it is found that the quantum conflnement gap 
severely interferes with alternative gap formation. Fig- 
ure ^ illustrates how the various gaps vary with stripe 
width, near xq = 0.25. The data are plotted vs aver- 
age doping, assuming a regular stripe array with 2-Cu 
wide AFM stripes (which do not contribute to the gaps) 
and N-Cu wide charged stripes of doping xo, giving an 
average hole doping x = Nxq/{N + 2). 

The CDW gap is most sensitive to stripe width, but 
in the narrowest stripes superconductivity is also sup- 
pressed (the suppression is stronger for a d-wave gap). 
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Strong instabilities are possible on the stripes, but they 
are shifted in doping away from xq — 0.25. Thus, the 
CDW instability requires both and Ej^^^ to be near 
the Fermi level; for a two-cell wide stripe, this is only pos- 
sible near x = 0. On the other hand, superconductivity is 
possible anywhere, if the coupling is strong enough. For a 
two-cell wide stripe, the optimal superconductivity arises 
when the Fermi level is at the one- dimensional VHS at 
the edge of one of the stripe subbands. This depends 
on t', and for f = -0.276i falls at x = 0.582 (a larger 
gap is found on the electron-doped side, x = —0.376). It 
should be noted that, even though the superconductivity 
is assumed to be d-wave, in general a finite minimum gap 
is found on the stripe, even when the CDW gap is zero. 
This is because the vanishing d-wave gap can be sampled 
only when the point {tt/2, tt/2) is sufficiently close to the 
Fermi surface, which in general requires N to be odd (re- 
call that the allowed values of kx are integer multiples of 
tt/{N + 1)) or very large. 
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+mV{\^\^+ < On >^ + <Tx>^ + <T >2), (16) 

where < T >^^< n >^ + < Ty >^ + < >^ and 
7fc = {cx — Cy)/2. The first two terms in Eq. [l^ renormal- 
ize the chemical potential and the hopping t respectively, 
and can be neglected. The terms in < > comprise a 
pseudospin triplet of CDW-like distortions, with rep- 
resenting a CDW similar to the one discussed above, Ty 
being related to the low-temperature tetragonal distor- 
tion, and Tz an orbital antiferromagnet, closely related 
to the flux phase. The remaining terms are a d-wave 
superconductor. The coefficients of the terms must be 
found self-consistently by solving the gap equations: 

^ < cl,C^,a >^1-X + 2{-ir < Tx >, (17) 

(7 

Im < cl^Ci+x,a >= -Im < 4^^Ci+y^a- >= (-l)*^' < Tz >, 

(18) 

Re < cl^c^+x,a >=< On> + <Ty> 

Re < c|_^c,+y,<, >=< On> - <Ty>, (19) 



FIG. 6. Gaps on a paramagnetic stripe, as a function of 
doping (equivalently: stripe width), assuming xq = 0.25, 
XcDW ~ 0.5eV, \a = 0.25eV. Squares — s-wave supercon- 
ducting gap; diamonds = CDW gap; triangles = combined 
gap; circles — total gap at (tt, 0), including the quantum con- 
finement gap. 



< clA+id >= ^- (20) 

The terms < 0„ > and < Ty > have pcecently been dis- 
cussed by Valenzuela and VozmedianoE3. A detailed dis- 
cussion of the competition between the three CDW-like 
modes is given in Ref. |6^. 



2. Electron-Electron Coupling 

In the above calculations, the A's arose from electron- 
phonon couplings. Similar contributions follow from 
electron-electron coupling, in an extended Hubbard 
model. For instance, the near-neighbor Coulomb repul- 
sion has the following mean-field expansion: 

V =4V^E4,.Cfe,- 
-2y <0„ >^(c. + c,)ct^c,,^ 

k^(T 



3. d-wave Superconductivity 

Retaining only the superconducting term in Eq. [l^ , 
the interaction can be derived from a quartic term 

with Vgj-= 2V{cos{kx — lx)a-\-cos{ky — ly)a). Assuming 
= IS.X cos kxa -\- Aj^ cos fcyO, the gap equations can be 
written in the form 

A, = ^ A,^,A, (22) 
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(i = x,y), with 

t&nh Ef/2kBT , ^ 
Aij = —2V 2_^cosliacoslja ^ , (23) 

r ' 

and E^- ^ ^{e^— ep)^ + Aj. 

For the uniform charged state (infinitely wide stripe) 
■Ax,x — ^y,yj and the gap symmetry can be simply an- 
alyzed. The symmetry can be either d-wave (Ay — 
—Ar^) or extended s-wave (Aj^ = +A3;), with the choice 
Ax^yAxAy > giving the largest gap. The A^^x term 
is always BCS-like, having the opposite sign from V, 
while the Aj. y term has the opposite sign from A^^^, 
the integral being dominated by the regions near the 
VHS's. Hence, there are two possibilities: (i) attrac- 
tive (y < 0) d-wave superconductivity or (ii) repulsive 
{V > 0) extended s-wave. However, the latter would re- 
quire \A.j;,y \ > \Ax,x\, which does not arise in the present 
model, so only case (i) is possible. These considera- 
tions readily generalize to a finite-width stripe, for which 

-^x^x ~r y,y rn 

Agterberg, et al.Ell recently introduced a model for 'ex- 
otic' superconductivity in multiband superconductors. If 
the Fermi surface consists of several inequivalent but de- 
generate pockets, the order parameter can consist of sym- 
metry allowed superpositions of the order parameters of 
the individual pockets. Equation ^ can be thought of as 
a form of exotic superconductivity, with the degenerate 
VHS's playing the role of hole pockets. 

B. Modifications due to Magnetic Order 




FIG. 7. Dispersion of linear antiferromagnetic (LAF) array 
along the Linear direction (X) (a), along with modifications 
due to d-wave superconductivity (b) or CDW order (c). U/t 
= 6, t'/t = 0, V/t = 2 (b), 0.1 (c). 

In the above calculations, it was implicitly assumed 
that the doping is high enough that the only role of the 



on-site repulsion U is to-jsenormalize the band parame- 
ters. However BaskaranE2 recently estimated that near 
optimal doping correlation effects remain stronger than 
the kinetic energy associated with hopping. Hence, it is 
important to look for stripe ground states which mini- 
mize this on-site repulsion (Class B stripes). The lin- 
ear antiferromagnet (LAF) stripes discussed in Section 
HI are a good candidate for the cuprate charged-atripes: 
they closely resemble the White-Scalapino stripesE^, have 
an appropriate doping, close to Xq = 0.25, include strong 
correlations, and have the additional advantage that a 
two-cell wide LAF charge stripe acts as a natural APB 
for the AF stripes. In this Section, we will explore these 
stripes, and show that they can be further stabilized by 
additional interactions. 

A special form of strongly correlated CDW is found to 
exist on a LAF. The charge and spin distribution is shown 
in the insert to Fig. ^, with the corresponding dispersion 
in Figure There is a strong antiferromagnetic ordering 
on one sublattice, while most of the holes are confined on 
the other, nonmagnetic sublattice. Whereas in a conven- 
tional CDW the charge density is zero on one sublattice 
and two on the other, in this strong coupling case the 
hole density varies from to 1, and there is no dou- 
ble occupancy, Fig. ||. Whereas the paramagnetic stripes 
were extremely sensitive to quantum confinement, these 
magnetic charged stripes are much less so: this CDW is 
stable almost independently of the stripe width. From 
Fig. it can be seen that the gapped Fermi surface still 
has hole pockets near (7r/2,7r/2), which would leadjio 
conducting stripes, consistent with optical propertiesEJ. 
However, it is only found near a hole doping x = 0.5, and 
so does not appear to be relevant for stripe physics in the 
cuprates. 
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FIG. 8. Linear antiferromagnetic (LAF) array with CDW, 
showing spin and doping distribution on different sites as a 
function of interaction strength V, with U/t — 6, t'/t — -0.276. 
Inset shows arrangement of atoms. 

Away from this doping, CDW instabilities are rela- 
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tively weak and it is possible to stabilize d-wave super- 
conductivity, Figs. |7|3,§. While the overall dispersion 
varies with stripe width, the superconducting gap is also 
relatively insensitive to the width, and actually increases 
for the narrowest stripes. Fig. Note that the order 
parameter is not a pure d-wave, the gap along the stripe 
being larger. Such a large anisotropy is not consistent 
with tunneling measurements of the gap; it is possible 
that the anisotropy is reduced by strong interstripe cou- 
pling. On the other tend, a large gap anisotropy has 
been found in YBCOEji-, where the stripes are aligned 
along the chain directionu. 




FIG. 9. Dispersion of a LAF with d-wave superconducting 
order for a uniform system (a) or a single stripe of width N = 
10 (b), 6 (c), or 2 (d) atoms. Darkness of line reflects relative 
intensity of dispersion feature. 
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FIG. 10. Linear antiferromagnetic (LAF) array with 
'd-wave' superconductivity, showing magnitude of gap along 
(j/) or across (x) the stripes, as a function of stripe width. 

Note in Figs. 0b, ^ that the combination of LAF and 
d-wave order leads to a finite minimum gap over the full 
Fermi surface. While the pure LAF phase is not stabi- 



lized by the VHS, the d-wave superconductivity is op- 
timized when the Fermi level is at the (tt, 0) VHS - at 
essentially the same doping, xq — 0.245, as the VHS on 
a paramagnetic stripe! 
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FIG. 11. Constant energy cuts of photoemission dispersion 
for a (2,2) stripe array, within 200 meV of the Fermi level. 
Lines = Fermi surface of bulk (or very wide) charged stripes. 
Relative intensity increases with darker shading, (a) Repre- 
sentative of single domain sample; (b) for multidomain sam- 
ple (symmetrized about the zone diagonal; (c) with diago- 
nal-suppressing matrix element, M = \cx — Cy\. 



VI. EXTENSION TO ARRAYS 
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FIG. 12. Constant energy cuts of photoemission dispersion 
for a (2,2) stripe array, with LAF charged stripes, within (a) 
30 (b) 100, (c) 200, or (d) 500 meV of the Fermi level. 

In Ref. the Fermi surface was calculated for a series 
of ordered stripe arrays. These results can now be com- 
pared to experimental photoemission dataLil. For this 
purpose we replot the data as integrated spectral weight 
over a finite energy cut within energy AE of the Fermi 
surface. Figure Elk shows a cut with AE = 200meV, ft 
the model of a 1/8 doped stripe array (i.e., x = 0.125)1 
The pattern is readily understood: the stripe superlattice 
leads to a number of quasi-one-dimensional bands; how- 
ever, due to structure factor effects they have significant 
intensity only near the original Fermi surface, solid line 
in Fig. [ll]a. For comparison with experiment, the cal- 
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culated spectral weight is symmetrized (0,7r) ^ (''■,0) in 
Fig. ^l|b to represent a sample with regions of stripes 
running along both X and Y directions. Finally, an 
empirical matrix element is included, Fig. [ll|c, which 
extinguishes spectral weight along the zone diagonal, 
(0,0) — ^ (tt, tt), similar to the matrix element assumed 
in analyzing Bi2212, Refs. The resuhing Fermi 

surface maps for several values of AE' are illustrated in 
Ref. [l^ for paramagnetic charged stripes, and in Figs. |l^ 
{x = 1/8) and 111 (a; = 0.1875) for LAF stripes. 
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FIG. 13. Constant energy cuts of photoemission dispersion 
for a (2,6) stripe array, with LAF charged stripes, within (a) 
30 (b) 100, (c) 200, or (d) 500 meV of the Fermi level. 

For both models, the stripe band nearest (tt, 0) is in 
good agreement with experiment: there is little disper- 
sion perpendicular to the stripe, while the intensity falls 
off toward (0,0) due to the structure factor effect. In 
general, the LAF stripes are in better agreement with 
experiment, since the additional subbands predicted for 
paramagnetic stripes (moving from (7r,0) toward (0,7r)) 
are not seen in the experiment. While the matrix element 
improves the agreement, theory suggests that this effect 
is present only for certain photon polarizationsO. One 
disagreement with experiment for both models is that for 
shallow energy cuts (30, 100 meV) the experiment still 
finds a smeared dispersion rather than a sharp Fermi sur- 
face. This is presumably an effect of stripe fluctuations. 

It should be noted that all the spectral weight in 
Figs. pT[]l^ is associated with the charged stripes; the 
lower Hubbard band of the AFM stripes lies below 0.5eV 
in LSCO. It is somewhat surprising that the spectral 
weight nearest the Fermi level is near (7r,0), since this 
is where the pseudogap arises. Nevertheless, our calcu- 
lation reproduces both the (quantum confinement) pseu- 
dogap, Fig. H, and the spectral weight distribution. 



VII. DISCUSSION 



A. Additional Evidence for Phase Separation 

1. Termination of Stripes 

The two most common interpretations for the pseudo- 
gap are in terms of either precursor pairing or .a .pom- 
peting order parameter, which|_may be magneticESEj or 
charge-density wave (CDW)^'n. We will argue in the 
next subsection that the latter possibility is more likely. 
Such competing instabilities arise naturally in a phase 
separation model (e.g., Refs. |77|,|9|), and we have long ar- 
gued that jth.e pseudogap is a manifestation of local phase 
separationlla. Itpshould be noted that a pseudogap arises 
in the nickelatesO in conjunction with stripe fluctuations, 
and turns into a true gap at the charge ordering temper- 
ature. Therefore, the fact that the pseudogap closes in 
the overdoped regime in the cuprates strongly suggests 
that the stripe phase terminates at the same^^inping. Di- 
rect evidence for this has recently appearedEnj'Ej'Ll. Fig- 
ure compares the intensity / of the inelastic neutros. 
peaks neaij-||(|7r, tt) as a function of doping in YBCOEil 
and LSCOHu. In a stripe picture, / should be a mea- 
sure of the fraction of material in AFM stripes. Re- 
markably, the intensity extrapolates to zero at nearly 
the same doping in both materials, even though Tc{x) 
peaks at substantially different dopings. For both mate- 
rials, this is the doping at which the pseudogap closes. 
While in YBCO, these peaks, i^ve been interpreted as a 
Fermi surface nesting effectOS, this cannot explain the 
charge order fluctuations, and in LSCO the Tranquada 
dataa is associated with elastic magnetic peaks. Further- 
more, a gtfidy of the neutron diffraction pair distribution 
functionEil for LSCO finds evidence for charge fiuctua- 
tions, presumably associated with stripes. The excess 
fluctuations are maximal near x = 0.15, and terminate 
near x — 0.25. Strong reductions of thermal conductiv- 
ity associated with-atripe scattering also terminate at a 
comparable dopingE3, while a recent optical studjitJ finds 
evidence for a quantum critical point at a similar |doping, 
X - 0.22. Moreover, in Bi2212, Tokunaga, et alO have 
introduced a new crossover temperature TmK based on 
Cu NMR, below which AFM correlations develop; they 
find TmK near x — 0.26. 

A recent NQR study of ths-slowing of spin fluctua- 
tions in RE substituted LSCOEa finds that the effective 
spin stiffness pl^^ (or equivalently the effective exchange 
constant) scales to zero at a comparable doping; the in- 
verted triangles in Fig. |l^ show 2it p^^ f / 'i&QK . As might 
have been anticipated from Sections II, IV, the doping de- 
pendence of p1^^ changes radically below x = 0.12. Note 
that while the integrated neutron intensity scales approx- 
imately with the area fraction of charged stripes, 2-Kpl^^ 
scales to ~460K as a; ^ 0. This is only 1/4 of the actual 
spin stiffness, 2-Kpg = 1.13J = ITBOif in the undoped 
AFM. The change by nearly a factor of four is suggestive 
of a dimensional reduction (lower coordination), but for 
an isolated straight spin ladder, a factor of two might 
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have been expected. 



2. Doping Dependence of Charged Stripes and 
Superconductivity 



* . 

0.25 J 



FIG. 14. (a) Maonetic inelastic scattering intensity "v^x 
in YBCO (squaresEJ) and LSCO (diamondsi and triangletl). 
Inverted triangles — (scaled) effective exchange constant 
in RE substituted LSCO, estimated from slowing of spin 
fluctuajtipnfS. (b) Corresponding Tc{x): solid line = YBCO 
(circlef£j); dashed line = LSCO. 

It is important to note the proximity of this termina- 
tion of the phase separation regime to the VHS: the ar- 
row in Fig. [l^ a showa-Jhe doping at which the pseudogap 
in the heat capacityEj closes, leaving an approximately 
logarithmic peakEj, while the x indicates the point at 
whick, photoemission finds the VHS crossing the Fermi 
leveEj. Termination of the stripe phase close to a VHS 
is an important prediction of our EPS model of stripes. 

Figure |lj has a bearing on the discussion given earlier 
of optimal doping in YBCO, or equivalently of the doping 
X at which the stripe phase terminates. Setting optimal 
doping at a; = 0.16 for all cuprates puts the termination 
at doping xq = 0.19. However, this value is clearly too 
low for LSCO: when the stripes are gone, magnetic cor- 
relations should be weak, but thetje is evidence for long 
range magnetic order at x = 0.2h and 0.21, the latter 
coupled with a suppression of TcCJ. On the other hand, 
both facts are compatible with xq ~ 0.25, Fig. |lj. 

The termination of the pseudogap has also been ii^tet- 
preted in terms of a quantum critical point (QCP)oeI. 
At first glance, it looks straightforward to distinguish 
a model of EPS, which is by nature first order, from 
a QCP, which is typically second order. However, the 
phases are actually quite similar. Due to strong coulomb 
effects, the EPS is restricted to nanoscopic scale, so that, 
e.g., the superconducting critical temperature Tc changes 
smoothly with doping. On the other hand, near a QCP 
with quenched disorder, long-timje-jCorrelations can lead 
to effective spatial inhomogeneitytJ. 



We suggested earlierB that the peak and hump fea- 
tures seen in photoemission from Bi2212 were associated 
with the charged and the AFM stripes respectively. As 
such, the intensity of the peak should have the dop- 
ing dependence predicted for charged stripes, with the 
intensity increasing from zero at half filling, approxi- 
mately linearlyj-with doping x. This has now been verified 
experimentallyOO. Moreover, the maximum intensity 
of the spectral weight occurs at the same doping xq 
discussed above, where the stripe phase terminates. Re- 
markably, the peak spectral weight closely tracks Tc, sug- 
gesting that the superconducting pairs U vq , cm the charged 
stripes, as predicted by several modela^Ea. Consistent 
with this, a number of measures of the strength of su- 
perconductivity (condensation oarx^iL critical current) 
are optimized at this same pointE3oE3 where the charge 
stripe intensity is maximum and AFM stripes vanish. 
The fact that Tc itself is actually optimized at a slightly 
lower doping may be a hint that stripes can enhance the 
superconducting gap, as found above. Fig. ^ 

The sharp falloff of spectral weight at higher doping 
may be an indication for a second regime of phase sepa- 
ration in overdoped samples, where the new phase is not 
intrinsically superconducting. Some additional evidence 
for such a phase separation is discussed below. 



3. Second Regime of Phase Separation 

It is important to note that as soon as the pseudogap 
closes, experiments find a peak in the density of states 
consistent with that expected for a VHS. This is discussed 
m RefJI on p. 1203 (Fig. 21); more recent evidence is in 
Ref. For even higher dopings there are hints of a new 
regime of phase separation. This is an improtant predic- 
tion of the model: if stripes are stabilized by lowering the 
electronic energy at a VHS, then there could be a second 
range of phase separation in the overdoped regime, with 
the non-superconducting phase totally unrelated to the 
AFM state at half filling. 

Earlier experimental evidence for this second regime 
of phase separation in the cuprates has been presented 
irL_B,eferences ^ (Section 11.6) and ^ (see also Wen, et 
alE3). Recently, LorarrO has presented a detailed anal- 
ysis of the heat capacity of La2_2:Sr2,Cu04 over an ex- 
tended doping range. In the range 0.26 < a; < 0.3, the 
temperature dependence of the Somerfeld constant 7 is 
consistent with that expected for a uniform phase close 
to a VHS (Ref. Fig. 21). For lower doping, 7 falls 
off too rapidly to be associated with rigid-band filling 
away from a VHS, but can be easily explained in terms 
of the opening oLa pseudogap. For higher hole doping. 



7 again falls ofitHEl too rapidly to be due to rigid band 



filling, although without an obvious pseudogap form- 
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ing. The near symmetry of 7 between the overdoped 
and underdoped sides is suggestive of two regimes of 
phase separation. Further possible evidence-for this sec- 
ond regime comes from tunneling studiesE£ll which find 
a subdominant order parameter in optimally and over- 
doped YBCO (i.e., in the same doping-aa,nge where a 
split superconducting transition is foundll23). T 
arise from an AUender-Bray-Bardeen mechanism 



lis. 
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could 
, with 

the d-wave superconducting stripes inducing s-wave (or 
dxy) superconductivity on adjacent norpial-Fermi-metal 
stripes. Also, microwave measurementsB find a crossover 
near optimal doping from superconductor-insulator to 
supercondjiietor-metal domains. 

UemurallHil has also recently proposed phase separation 
in the overdoped cuprates. In extending the Uemura plot 
to the overdoped regime, a 'boomeraugl effect is found: 
both Tc and n^/m decrease in paralleL£j. Here Ug is the 
superfluid density and m is the effective mass, and their 
ratio is extracted frcpa-JiSR measurements of the penetra- 
tion depth. Uemurall23 suggests that this proportionality 
between and jig/m can most easily be accounted for in 
terms of (nanoscale) phase separation between the super- 
conductor and a non-superconducting pliaae. Thus for 
both overdoped and underdoped cupratesE2j, the scaling 
between Tc and can be understood in terms of phase 
separation. This has important implications for super- 
conductivity: superconductivity arises only at a unique 
composition near optimal doping. The present results 
suggest that even at this doping, the superconductivity 
must compete with a second order parameter. 



B. Superconducting Fluctuations 

Above, we have presented evidence that the optimum 
Tc occurs at a different doping in LSCO than for the other 
cuprates, while the stripes terminate at approximately 
the same doping for all cuprates. Further evidence for 
this can bo-Jo und I py comparing the superconducting 
fluctuationalHZl^llijO, Fig. 15, In all three materials. 



LSCO, YBCO, and Bi2212, some superconducting fluctu- 
ations appear to persist in the low doping regime, to tem- 
peratures in excess of optimal Tc (although considerably 
less than the pseudogap temperature). The fluctuating 
superconductivity in LSCO is very similar to that found 
in the other cuprates. On the other hand, if one assumed 
a universal curve of superconductivity, with Tc/Tc^max 
vs X the same for all cuprates, then the scaled supercon- 
ducting fluctuations in LSCO would be at anomalously 
high temperatures. It seems simpler to assume that in 
LSCO long range superconducting order is_aiiomalously 
suppressed (as also suggested by Baskarantill) , in which 



case there is no reason for optimal Tc to fall at the same 
doping. 

In Bi2212 and YBCO, the onset of superconduct- 
ing fluctuations falls close to the 'strong pseudogap' 
regimdliSEHl, and hence at temperatures substantially 



below the weak pseudogap temper aturellijili3, making it 
unlikely that the weak pseudogap is associated with pre- 
formed pairs. In LSCO, there is some confusion, since the. 
pseudogap temperatures reported by Batlogg, et aLEis 
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are considerably higher than those found by Ido, et al 
The Ido r*'s fall close to the fluctuation onset tempera- 
tures in Fig. nB. 




0.2 

X 



FIG. 15. Doping dependence of Tc (solid line for YBCO, 
Bi2212, dashed line for LS CO) and superconducting fluctua- 

11, LSCO 



110 



and (diamonds}^ 



tions in Bi2212 (circles) [ |1071] 
(triangles) f LOS], and YBCO (closed squares) [ lOE] and (open 
squares) 



The fluctuations observed in Fig. nSl are very sugges- 
tive of the theoretical results of Fig. |10|: thus, quantum 
size effects enhance the average gap on a narrow stripe, 
which will lead to an enhanced mean fleld transition tem- 
perature in the underdoped regime. However, the cou- 
pling between stripes also weakens, leading to a loss of 
interstripe coherence, rSO, the pairing is only evident as 
enhanced fluctuations Jlii 



VIII. CONCLUSIONS 

Recent experiments have provided considerable evi- 
dence for the presence of stripes and EPS in the cuprates, 
but there remain many questions of how universal these 
are, how they arise and vary with doping, and how they 
interact with superconductivity. We have here elaborated 
our earlierB model of stripes driven by frustrated phase 
separation, in particular adducing evidence that the dop- 
ing on the charged stripes is close to x ~ xq ^ 0.25, 
and that when the average doping approaches this value 
EPS terminates. Moreover, near Xcr = xq/2 there 
is a crossover in stripe properties: for x < Xcr the 
charged stripes are quantum confined, for x > Xcr the 
AFM stripes are so confined. This model can explain 
the 1/8 apomaly {x = Xcr): the anomalous Hall effect 
{Rh 0]B for a; < 1/8 (charged stripes confined, hence 
one-diiriensional) , and the growing spin gap in YBCO for 

X > i/m. 

On the important issue of the structure of a charged 
stripe, we have explored a number of possibilities without 
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coming to any final conclusions. While there is evidence 
that superconductivity lives on the charged stripes, there 
also appears to be a second instability on these stripes, 
which stabilizes the stripe phase while competing with 
superconductivity. Wc have shown that a semiquanti- 
tative understanding can be achieved by looking at the 
properties of a single doped ladder, and we have discussed 
how a number of instabilities (both CDW and cuper- 
conducting) vary with ladder width. We showed that 
strong correlation effects could lead to charged stripes 
with a residual magnetic order, introduced a simple 
model for White-Scalapino-likc stripes, and found novel 
superconducting and CDW instabilities associated with 
such stripes. We illustrated how stripe order would af- 
fect ARPES spectra, both dispersions and Fermi surface 
maps. Future studies will apply the model to describing 
other properties of the cuprates. 

Certain anomalous features of strong-correlation cal- 
culations may find an explanation in underlying phase 
separation. Thus, the vanishing of the renormalized hop- 
ping parameter t near half filling in slave boson calcula- 
tions may reflect the vanishing of the charged stripes at 
half fiUingEl, while the frequently-observed pinning of the 
VHS near the Fermi levelBEill is consistent with VHS- 
stabilized charged stripes. 

Finally, the idea of a commensurability effect near 1/8 
doping, leading to a coexistence of domains for x > 1/8, 
provides a simple explanation for a large variety of ex- 
perimental findings, including the saturation of the Ya- 
mada plot, the direct observation of domains in STM 
studies, and a variety of microwave anomalies. This may 
also lead to a resolution of the combined puzzle of mag- 
netic neutron scattering incommensurability and the neu- 
tron resonance peak. A stripe model provides a natu- 
ral explanation of the incommensurability for x < 1/8, 
including a stripe reorientation transition at the metal- 
superconductar transition near x 0.053. However, a 
band pictureEJ (with Ep close to a VHS) provides a su- 
perior model for the combined, frequency-dependent in- 
commensurability cum resonance peak found near opti- 
mal doping in YBCO and Bi2212. An EPS crossover from 
stripes to domains near a; ~ 1/8 would provide a natu- 
ral explanation pof these phenomena. Interestingly, the 
only calculationEa to approximately describe the doping 
dependence is based on a slave boson model, which as 
noted above tends to mimic EPS behavior as a; — > 0. 
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